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Abstract We present new optical spectroscopy of the 7-ray binary LS 5039. Our 
data show evidence for sub-orbital modulation in the radial velocities with ampli- 
tude ^7 km/s and period ^ Porb/4. This short-term oscillation is stable over at least 
7 years and it is likely triggered by non-radial oscillations of the 06.5V optical star 
We also present the results of a spectroscopic campaign on MWC 148, the optical 
counterpart of the new 7-ray binary candidate HESS J0632H-067. Long-term vari- 
ations in the Ha and Hp emission line parameters are clearly detected which, if 
modulated with the binary orbit, would imply a period >200 days. 



1 Introduction 

There are only three confirmed galactic High Mass X-ray Binaries (HMXBs) with 
persistent TeV emission (PSR B1259-63, LS 5039 and LSI +61 303, see [1]) and 
one recently proposed candidate (HESS J0632H-067 [2]). The VHE emission in the 
3 confirmed binaries is strongly modulated with the orbital period, suggesting that 
the emitter is rather compact and close to the massive star. However, the origin of 
the VHE emission remains unclear, with competing leptonic and hadronic scenarios 
which invoke either inverse Compton scattering of stellar photons or proton-proton 
collisions in a jet/pulsar wind scenario. Obviously, the nature of the compact star 
is an important ingredient in the different 7-ray production models. PSR B 1259-63 
contains a pulsar whereas a black hole is not ruled out in LS 5039 and LSI +61 
303. In this context optical studies can bring new insights through constraining the 
compact object mass with dynamical studies. 
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2 Revised Orbital Solution in LS 5039 

LS 5039 was the first 7-ray binary detected by HESS, with a TeV luminosity of 
-10^"^ erg s-K It was previously classified as a HMXB with spatially resolved ra- 
dio emission at milliarcsec scale, atributed to either a relativistic jet or a pulsar 
wind. The first determination of the 3.9d orbital period and system parameters was 
reported by [3]. The orbital parameters were found to depend on the spectral lines 
used in the cross-correlation analysis. In particular, the Balmer and Hel solutions are 
blueshifted with respect to the Hell solution and also show a larger scatter This is 
expected if low excitation lines are contaminated by P-Cyg profiles and, hence, the 
Hell velocities provide a much better description of the orbit of the optical star [3]. 
A tentative low inclination angle in the range 23-27° was derived under the assump- 
tion that the companion star is pseudo-synchronized, i.e., synchronized at periastron. 
This results in a compact object of 3-5 M© and, thus, the possibility that LS 5039 
harbours a low-mass black hole (see [3] for details). 
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Fig. 1 Updated radial velocity curve in LS 5039. Left: velocities folded on the best ephemeris with 
a two Fourier term fit sumperimposed. Filled circles indicate INT data from 2002 & 2003 reported 
in [3], open triangles are SAAO data from 2005 and open circles INT data from 2007 and 2009. 
Right: Same as left but with an extra sine wave to account for the short-term modulation. 

New spectrocopy was obtained in 2005, 2007 and 2009 using the 1.9m telescope 
at SAAO and the INT at La Palma. Fig. 1 displays the new radial velocity points, 
together with data from [3], folded on the best orbital solution. The orbital modula- 
tion is clearly distorted by a short-term oscillation with amplitude ^ 7 km s^' and a 
~ld period or Porb/4. In order to test the significance of this perturbation we fitted 
the eccentric orbital modulation with a 2-term Fourier series and obtain a reduced 

of 11.3 (left panel in Fig. 1). We subsequently added an extra sine wave to ac- 
count for the short period oscillation and find that the reduced drops to 7. 1 (right 
panel in Fig. 1), thus confirming that the short-term modulation is highly significant. 
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We note that this modulation must be stable over our 7 years database, otherwise it 
would be blurred in the phase folded radial velocity curve. 



Table 1 Orbital solutions for LS 5039 



Parameter 


Eccentric Fit 


Eccentric Fit 






+ Id oscillation" 


Poib (days) 


3.90597± 0.00009 


3.90608± 0.00008 


To (HJD-2450000) 


2478. 11± 0.08 


2478.08± 0.06 


e 


0.30± 0.03 


0.35± 0.03 


wO 


215± 6 


212±5 


7 (km s^') 


17.3± 0.5 


17.3±0.5 


Kx (kms-') 


23.1±0.9 


24.2± 0.9 


a\ sin! (R0) 


1.70± 0.07 


1.75± 0.07 


/(M) (Ms) 


0.0043± 0.0005 


0.0047± 0.0006 


rms of fit (km s^ ' ) 


8.3 


6.4 



° The period has been fixed to Poit^^ 



Ignoring this perturbation can certaintly bias the final orbital solution as illus- 
trated in Table 1 . The 2 columns compare the binary parameters obtained through a 
proper eccentric model fit to the Hell velocities with and without the Id oscillation. 
Both, the eccentricity and radial velocity semi-amplitude K\ are underestimated if 
the Id oscillation is neglected. Also the periastron angle and the absolute phasing 
differ The key binary parameters are still consistent within \-a but the difference 
can be larger when the fit is performed over scarcely sampled radial velocity curves. 
We note that this short-term modullation is very reminiscent of the 5 km s^' oscil- 
lation detected in Vela X-1 at 1/4 of the orbital period [4]. Following the work on 
Vela X-1 we propose that the Id modulation in LS 5039 is caused by tidally excited 
non-radial oscillations of the 06.5V star in its eccentric orbit. The presence of this 
oscillation may also have an impact on the VHE radiation models because of the 
role played by the stellar photons and hence it should be investigated. 



3 Probing Binarity in MWC 148 

MWC 148 is the optical counterpart of the point-hke TeV source HESS J0632H-067 
and it has been recently proposed as a candidate 7-ray binary based on the spectral 
properties and variability observed in both radio and X-rays [2]. In Oct 2008 we 
started a spectroscopic campaign to observe MWC 148 using several telescopes at 
the Teide and Roque de los Muchachos observatories. The averaged spectrum shows 
Ha and emission lines superimposed on a BO-type stellar spectrum dominated 
by Hel absorption lines (without any emission component). From the Hel profiles 
we estimate a rotational broadening vsin/ = 370 km/s. Since Be stars are known 
to rotate at typically 0.7-0.8 times the critical rotation velocity (i.e. ^565 km/s 
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for a BOV star) the binary inclination must be moderately high, in the range ^55- 
70°. Cross-correlation of the absorption lines with a broadened template star gives 
velocity excursions of a few tens km s^' but no clear periodicity is observed. 



Fig. 2 Time evolution of 
emission Une parameters in 
MWC 148 from October 2008 
until April 2010. Top panel 
shows the EW, middle panel 
the FWHM and bottom panel 
the centroid velocity obtained 
through a Gaussian fit. The 
same trend is observed in Hj3 
but with lower amplitude. 




so GO 



However, a dramatic change in the emission line parameters is detected between 
November 2009 and January 2010 (see Fig. 2). The lines become narrower, the EW 
drops by a factor 2 and radial velocity variations are also seen. If these changes are 
modulated with the binary phase then the orbital period must be longer than 200 
days. However, superorbital Ha variations due to precession of the circumstellar 
disc are also commonly seen in Be X-ray binaries [5]. Clearly many more observa- 
tions are needed before the orbital period and binary parameters can be constrained. 
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